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Abstract 

The Karman-Howarth theorem is derived for the Lagrangian averaged Navier-Stokes alpha 
(LANS— a) model of turbulence. Thus, the LANS— a model's preservation of the fundamental 
transport structure of the Navier-Stokes equations also includes preservation of the transport 
relations for the velocity autocorrelation functions. This result implies that the alpha- filtering 
in the LANS— a model of turbulence does not suppress the intermittency of its solutions at 
separation distances large compared to alpha. 

1 Introduction 

Karman-Howarth (KH) theorem 

The dynamics of the two-point velocity autocorrelation functions in turbulence contains informa- 
tion complementary to the spectral description. These autocorrelation functions directly connect 
the concept of scale with the result of an actual flow measurement. However, the autocorrelation 
functions yield no information about the energy, or enstrophy, that is contained in their interval of 
separation. Instead, the third-order correlation functions do provide information about the fluxes 
of energy and enstrophy as functions of separation distance. Thus, the spectral and autocorrela- 
tion dynamics give complementary information about the same phenomenon. The KH theorem 
relates the time derivative of the two-point velocity autocorrelation functions to the divergences 
of the third-order correlation functions. Perhaps not surprisingly, because of its fundamental im- 
portance, the line of investigation that began in 1938 with the KH theorem is still being actively 
pursued. 

The invariant theory of isotropic turbulence was introduced by von Karman & Howarth (1938) 
and refined by Robertson (1940), who reviewed the Karman-Howarth (KH) theorem in the light 
of classical tensor invariant theory. The physical importance of the KH theorem in turbulence 
modeling is undeniable. According to Monin & Yaglom (1975) (vol II, page 122) the KH theorem's 
dynamical equation for the two-point autocorrelation function of the fluid velocity, "plays a basic 
part in all subsequent studies in the theory of isotropic turbulence." A homogeneous (but not 
necessarily isotropic) version of this theorem was introduced by Monin (1959). This non-isotropic 
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variant of the KH theorem was discussed further by Monin & Yaglom (1975) (vol II, page 403), 
but there was a gap in the proof. Correct proofs were given independently by Frisch (1995) and 
Lindborg (1996). These proofs are reviewed in Hill (1997), who concentrates on the logical steps 
needed to eliminate pressure- velocity correlations in the KH theorem without assuming isotropy. 

The isotropic Karman-Howarth theorem was extended to magnetohydrodynamics (MHD) in 
Chandrasekhar (1951) and was recently analyzed further using Elsasscr variables by Politano &i 
Pouquet (1998). Extensions of the KH theorem and its corollary - the "—4/5 law" of Kolmogorov 
(1941b) - to include passively advected scalars such heat and chemical concentration (or buoyancy) 
were given by Yaglom (1949) and were later refined by Antonia, Ould-Rouis, Anselmet Sz Zhu 
(1997), who also compared with Yaglom's results for the case of stratified fluids. The effect of 
fluid helicity on the KH theorem was discovered in Chkhctiani (1997) and analyzed in terms 
of structure functions in Gomez, Politano &: Pouquet (1998). Arad, L'vov & Procaccia (1997) 
extended the fundamental results in Robertson (1940) by considering projections of the fluid 
velocity autocorrelation dynamics onto other irreducible representations of the SO (3) symmetry 
group, besides the scalars under rotation. Modern computational resources are frequently applied, 
such as in Fukayama et al. (2000) and references therein, for numerically studying intermittency 
in turbulence, which is represented by the anomalous scaling behavior of its higher order structure 
functions, relative to those in the KH theorem. 

Navier-Stokes alpha model 

Recently, a modification of the Navier-Stokes (NS) equations known as the Navier-Stokes alpha 
model was introduced and its solution behavior was compared with experimental data for turbu- 
lent flow in pipes and channels in Chen et al. (1998, 1999a, 1999b). These modified equations 
essentially filter the fluid motion that occurs below a certain length scale (denoted as alpha), which 
is a parameter in the model. Amongst other differences from traditional turbulence modeling ap- 
proaches, this alpha-filtering approach differs from the large eddy simulation (LES) approach by 
preserving the basic transport structure of the exact NS equations. Direct numerical simulations 
(DNS) of this model for forced homogeneous turbulence were performed in Chen et al. (1999c). 
Homogeneous turbulence decay was also simulated numerically using this model in Mohseni et 
al. (2000). Both these simulation studies showed the NS— a model to be considerably more com- 
putable than the exact NS equations, while preserving essentially the same solution behavior as 
NS at length scales larger than alpha. A review of the basic properties of the Navier-Stokes alpha 
model and its early development is given in Foias, Holm and Titi (2001a). See also Foias, Holm 
and Titi (2001b) and Marsden and Shkoller (2001) for recent analytical and geometrical results 
for this model. 

The alpha model's original derivation (without viscous dissipation) in Holm, Marsden and 
Ratiu (1998a, 1998b) was motivated by seeking a generalization to n dimensions of the integrable 
one-dimensional shallow water equation due to Camassa and Holm (1993). This derivation was 
based on the Lie algebraic and geometric properties of Hamilton's principle for fluid dynamics in 
the Eulerian description. However, motivated by its later derivation as a Lagrangian averaged 
fluid model in Holm (1999), Holm (2001), and by a related derivation in Marsden and Shkoller 
(2001), we shall henceforth refer to this model as the Lagrangian-averaged Navier-Stokes alpha 
model, or the LANS— a model. 

Although several analytical and geometrical results are known for the LANS— a model, many 
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open questions still remain. An outstanding question concerns the effect of alpha on intermit- 
tency in turbulence. One supposes that filtering at scales smaller than alpha would suppress 
intermittency on those scales. However, what is the effect of alpha-filtering at small scales on 
intermittency in the larger scales? Proving the KH theorem for the LANS— a model is a required 
first step toward addressing this outstanding question. Remarkably, the effects of the alpha- 
filtering are confined to separations of order 0{a) or smaller, in the analogs of the KH theorem 
and in Kolmogorov's energy balance law proved here for the LANS— a model. 

The KH theorem for the LANS— a model proved here provides an exact result that demon- 
strates how the introduction of the length scale alpha affects the dynamics of the velocity auto- 
correlations as a function of separation between two points fixed in the domain of an isotropic 
LANS— a flow. These effects turn out to be negligible at separations that are large compared to 
the flltering scale alpha (r ^ a). To the extent that the same result persists in the higher correla- 
tion functions and structure functions for the LANS— a model, one may expect the intermittency 
of the alpha model to be the same as that for Navier-Stokes at large separations (r ^ q). 

The LANS— a model equations 

The motion equation for the unforced LANS— a model is given by 

+ u • Vv + UjVn^' + V(^p- ^|up - y |Vup^ = uAv , 
where v = u — a^Au , V • u = . 

The left side of this motion equation preserves the transport structure of ideal fluids, especially 
the Kelvin-Noether theorem, see Holm, Marsden and Ratiu (1998a, 1998b). Here, the specific 
momentum v is related to the fiuid velocity u, via the Helmholtz operator 1 — A which contains 
the length scale a. Although derived from a different approach, these equations are formally 
similar to the equations for the motion of a 2"'^ grade fluid, but with a different dissipation 
operator. Following Camassa & Holm (1993) for the original one-dimensional version of this 
equation, the motion equation for the alpha model may also be expressed equivalently in a way 
that emphasizes its nonlocality, as 



(l-a^A)(^atU + u- Vu-i/Auj +Vp = -divr, (1.1) 
where r = Vu • Vu + Vu • Vu^ - Vu^ • Vu , V • u = . 

Similar forms of the motion equation for a 2"°' grade fluid also appear in Dunn & Fosdick (1974), 
modulo a crucial difference in dissipation.]^ Upon inverting the Helmholtz operator in three 

^The formal similarity of the LANS— a model with the 2"'' grade fluid equation (8.1) on page 241 of Dunn & 
Fosdick (1974) resonates with a remark made already in von Karman & Ifowarth (1938) page 197, "The expression 
. . . for the correlation tensor is exactly of the same form as that for the stress tensor for a continuous medium 
when there is spherical symmetry." Mathematically, the crucial difference betw een the LANS— a, and 2"^ grade 



fluid models lies in their different dissipation operators. In the motion equation (1.1) for the fluid velocity, u, the 
LANS— a model has diffusive viscous dissipation — z/Au, introduced in Chen et al. (1998, 1999a, 1999b) for their 
com parisons with turbulent flow data, while the 2"'' grade fluid has only damping. That is, — vAu in equation 



(1.1) is replaced by fiu in the motion equation for 2"'' grade fluids, for a positive constant /i. Physically, there 



are also significant differences in the interpretations of both the solutions and the parameters that appear in the 
equations for the two models. 
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dimensional Cartesian components, this becomes 

dm + dk ( Uiu'' +p6f)- i^Aui = -a^ dk {ga * rf^ , where i,j,k £l, 2, 3, 

Here (1 — a^A)p = p and ga is the Green's function for the Helmholtz operator in the flow domain 
with boundary conditions compatible with satisfying NS on the boundary. The free space Green's 
function for the Helmholtz operator in three dimensions is gaii") = exp — (r/a), which is the 
well-known Yukawa potential for r = |x' — x|. This free space Green's function should apply to 
the extent that the turbulence is isotropic (so that a can be taken as constant) and away from 
boundaries at distances greater than order 0{a). 

This form of the LANS— a equation demonstrates that the nonlinear stresses proportional 
to are filtered by the Green's function gair) of the Helmholtz operator, whose filter width 
is of order 0{a). Given a, the term g^ * rf applies at separations r such that a/r < 1 and is 
negligible for a/r <C 1. For those larger separations, the LANS— a equation reverts to NS. The 
alpha-filtering is a regularization. With it, the solutions of LANS— a are well-posed and have a 
finite dimensional global attractor under L2— bounded forcing, as shown in Foias, Holm and Titi 
(2001b). See also Marsden and Shkoller (2001) for a proof of their well-posedness in the presence 
of boundaries. Such well-posedness results are not known to hold for the NS equations. However, 
solutions of the LANS— a equations do converge to solutions of the NS equations uniformly as 
a — > for any positive viscosity, as also shown in Foias, Holm and Titi (2001b). 

The equations governing the correlations 

In deriving the correlation dynamics, one may regard {ga*T'j ) as a subgrid scale (actually, sub— a 
scale) stress tensor arising from the alpha-filtering procedure represented by convolution with the 
Green's function ga{f), which decreases exponentially in separation with scale length a. See Holm 
(1999), Holm (2001) and Marsden and Shkoller (2001) for further analysis and explanation of how 
this viewpoint arises, upon applying Lagrangian averaging. We denote u{x',t) = u' and begin 
our investigation of the correlation dynamics by computing the ingredients of the partial time 
derivative dt{viUj), 

dm + dk ( Vivf' + p5^ - C? Ui^raV^^^ ) = vAvi , 

dtu' j^dl^u' ju'^ ^p' 8]) + a^dk (5*r')) =vA'u\. 

Wc cross multiply and add these equations, average the result ( • ) and use statistical homogeneity 
in the following form, with ^ = x' — x, in the traditional KH notation, 

— n=^n=-— n 

de ^ ' dx'^ ^ ' dx^^ 



to find 



( ) - ^ ( (^^^' - Ui,mU^^'^)u'j ) + ^{{vip')6^-{ u'^p) 6^ ) 

d 



( Vi {u'^u'^ + a2 5 * T'^^ ) ) = 2Z/A4 ( {viu' ^) 



D. D. Holm 



KH theorem for LANS— a model 



May 16, 2001 



5 



where is the Laplacian operator in the separation coordinate ^. Next, we symmetrize in i,j 
and use the relation obtained from homogeneity, 

( ViU' jU'^ + VjU'iU'^ ) = — [v'iUjU^ + v' jUiU^ ) , 

to find the homogeneous correlation dynamics for LANS— a, 

dt {viu'j + Vju\) - (t^j - a" 4 - nf^.) = 2v/:^^[v^u'j+Vju'i) . (1.2) 
In this equation, the symmetric tensors T^, H^j and S^j are defined as 

Tjj = ( ViU 'j + Vju'i + v 'iUj + V 'jUi ) u'' , 

5fj = ( Ui^mU 'j + Uj^rnU ' i)u'^'^ + [vi g * T '] + Vj Q * T ) . 

Imposing isotropy 

We now suppose the LANS— a solution is isotropic and follow the approach of Karman & Howarth 
(1938), as refined by Robertson (1940) and Chandrasekhar (1951) using the invariant theory of 
isotropic tensors. By a standard argument, isotropy implies we may drop the pressure-velocity 
tensor 11^^. Hence, we rewrite equation ( |1.2D as 

I 2., = {rt, - a' 4) + 2.A,Q,, , (1.3) 



with the corresponding definition, Qij = {viu'j + Vju'i). According to their definitions, the 
tensors Qij and Tjj are both symmetric and divergence-free in their indices i,j for constant 
a. For consistency with the isotropy assumption, equation ( |1.3D implies that Sfj must also be 
symmetric and divergence- free in its indices 

Remark. Thus, DNS could check whether a 7^ is consistent with isotropy in a given flow 
regime by checking whether the divergences of S^j vanish in its indices 

According to the classical theory of invariants discussed in Robertson (1940) and Chan- 
drasekhar (1951), these three symmetric, divergence-free, isotropic tensors may each be expressed 
in terms of a single defining function. In particular, the isotropic u — v autocorrelation tensor is 
given by 

Qij = curl (Q Eiji^^ ) = rQ' - 5,,) - 2Q% , 

with defining function Q{r, t) and Q' = dQ/dr in the KH notation. The isotropic triple correlation 
tensor is 

7:^^ = curl (r {iiSjM + ijE.u ) ) 

= \t 'iiijik -{rT' + 3T) (^Cidjk + ^jSik) + 2T6ijCk , 
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with defining function T(r, t) and antisymmetric tensor Eiji . Hence, we compute tlie divergence, 

^ Ti'; = curl ((rT' + 5T)e,,,e) = curl (l(r5r) 'e,,,^^ ) . 

This is formula (45) in Chandrasekhar (1951) and, of course, it agrees with the corresponding 
formula (4.13) in Robertson (1940). Likewise, the isotropic mean sub— a scale stress tensor S^j 
must also take the same form, 

^ S^j = curl ((r^' + 5S)eijeC' ) = curl (^(r^S) 'eijiC' ) , 

with defining function S{r,t). 

Remark. Because to the presence of curl in their definitions, the defining functions T and a'^S 
have dimensions of energy dissipation rate, {v?)'^/'^ /r. Their relations to the defining functions in 
Karman & Howarth (1938) will be discussed in a later section. 

The equations governing the defining scalars 

According to Robertson (1940), the scalar defining the Laplacian of a second order isotropic tensor 
is obtained by operating with 

V dr^ r dr J dr dr ' 
on the scalar defining the original tensor. That is, 

D(Q)=r-\r^Q')' . 

The scalars defining the various second order tensors in equation ( |1.3| ) are, therefore. 

As Robertson (1940) points out, upon assuming isotropy, the tensor equation (ll.3|) is entirely 
equivalent to the corresponding scalar equation. This observation proves the following: 

Theorem 1.1 (KH theorem for the LANS— a model) Let the LANS— a model flow be isotropic. 
Then, the scalar equation 

^ = + 5) (T - a^S) + 2uD{Q) , (1.4) 

defines an exact relation among the 2nd and 3rd correlation functions, and the stress tensor S^j. 

Remark. Formula ( |1.4| ) is the analog for the LANS— a fluid equations of the KH theorem for NS 
turbulence. 

Upon setting dQ/dt = —2^/3 for the energy dissipation rate in three dimensions and drop- 
ping the viscous terms in equation ( |1.4D - as is appropriate for separation r in the inertial range 
- one finds the energy balance relation for the LANS— a model. 
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Here, denotes the average dissipation rate of the total kinetic energy for the LANS— q model, 
given by Ea = ^/u • vd^x. Integration of the energy balance relation (1^) then proves the 
following: 

Corollary 1.2 ("—2/15 law" for the LANS— a model) In the inertial range and for arbi- 
trary ratio a/r, the LANS— a model satisfies the "~2/15 law," 

-^e^ = T-a^S. (1.6) 
15 

Remark. Formula ( |1.6D is the analog for the LANS— a fluid equations of Kolmogorov's "-4/5 
law" for NS turbulence. 

Comparisons with NS turbulence theory for a/r ^ 

Foias, Holm and Titi (2001b) show that solutions of the LANS— a model converge to solutions of 
the NS equations as a — > uniformly for any positive viscosity. Therefore, to compare the KH— a 
theorem with the results of Karman & Howarth (1938) for the NS equations, we may consider 
the limit as a/r 0. In this limit, one may neglect a'^S in equations ( p..5[ ) and ( |1.6D . We follow 
Robertson (1940) in identifying the KH double correlation scalars f{r,t), g{r,t) as 

(^)/ = Q, and {^)g = \rQ'. (1.7) 

Likewise, the KH triple correlation scalars h, k and q are identified in terms of T(r, t), as {u^ 
rT/2 for h, as well as. 



{u^fl^k = -\ r s^Tds and (u^f/^q = T s^T 
Jo Sr* _/o 

for k and q. Using the relation T = 2{u'^)^/'^h/r yields, 



r 

ds — - T . 
4 



Hence, when a/r — > 0, equation ( |1.4D of the KH— a theorem LI yields equation (51) of Karman 
& Howarth (1938), namely, the KH equation. 



dt \dr r 



2(n2)3/2/i + 2i/(n2)^ 
or 



[1.9) 



One may refer to Monin & Yaglom (1975) page 122, for the KH equation in their notation. When 
the factor 1/6 relating third order autocorrelation functions and structure functions is introduced, 
this is also equation (3) of Kolmogorov (1941b), leading to the Kolmogorov's —4/5 law for Navier- 
Stokes fluids. See also Landau & Lifschitz (1987) for additional discussion of this fundamental 
result. 



Thus, the limit a/r ^ of formula (L4) of the KH— a theorem recovers the expected classical 



results for homogeneous, isotropic, NS turbulence. 
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2 Conclusions 

Main results. The main results of the paper are, as fohows. 

• Equation ( |1.2| ) for the LANS— a dynamics of the velocity autocorrelation functions with 
homogeneous statistics shows that the alpha model's preservation of transport structure 
extends to preserving the form of these transport equations, as well, modulo a shift in the 
term involving third moments that accommodates the contribution from the sub-alpha scale 
stress tensor. 

• The modifications of the KH theorem and Kolmogorov's energy balance law for NS to include 
the effects of alpha-filtering in the LANS— a model are relatively simple and straightforward. 

• Equation ( |1.4| ) in the KH— a theorem gives the alpha-analog of the Karman-Howarth 
theorem for isotropic turbulence. 

• As a corollary of the KH— a theorem, equation ( |1.6| ) gives the "—2/15 law" alpha- 
analog of Kolmogorov's "—4/5 law" for energy balance in the inertial range. (The 
relative 1/6 in these coefficients arises from the 1/6 relation between autocorrelation 
functions and structure functions in isotropic turbulence.) 

• The second term in the —2/15 law in equation (|1.6| ), the 5— term which is proportional 
to a^, is very reminiscent of the quantity that appears in the corresponding —2 law 
for enstrophy cascade in 2D turbulence. The latter expression contains two powers of 
vorticity and one power of velocity. For example, see the Appendix B of Eyink (1996), 
where the identity is derived in detail. In the same way, the 5— term in (|1.4D contains 
two powers of velocity gradient and one of velocity. It should thus be the dominant 
term (compared to the first T— term) when r < a and the velocity is nearly smooth. It 
then corresponds to an "enstrophy- like" cascade, in agreement with the considerations 
of Foias, Holm and Titi (2001a) that determine the spectrum in that rangej^ 

These exact results confirm that alpha-filtering affects the homogeneous isotropic statistical 
properties of NS only at separations r of order 0(a), or smaller. This is because of the exponential 
decrease with separation in the Green's function for the Helmholtz operator. The same exponential 
decrease with separation applies for the dynamics of the higher order autocorrelation functions 
and structure functions, whose anomalous scaling properties for NS indicate the intermittency of 
turbulence. Consequently, one may expect the intermittency of the alpha model will be the same 
as for NS at large separations, for which r/a ^ 1. This conclusion may apply more generally for 
other types of filtering, as well. 

Some aspects of these results were not entirely unexpected. After all, one could enforce 
homogeneity and isotropy on the velocity correlation functions for any fluid model (e.g., an LES 
model, or the 2"'^ grade fluid model). One could then compute its corresponding KH equation 
by the same classical methods as those used here. In general, this procedure could introduce 
complications, such as requiring several more scalar deflning functions. Remarkably, the KH— a 
theorem for the case of constant alpha may be written at the cost of just one additional scalar 
function, S(r,t), arising from the sub-alpha scale stresses. The LANS— a model preserves the 
fundamental transport structure of the NS equations. It also preserves the form ( |1.3D of the 

am grateful to G. Eyink for this observation. 
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transport relations for the velocity autocorrelation functions. Thus, the interpretations of its 
consequences turned out to be relatively straightforward in this case. 

Helical turbulence. These results readily extend to helical turbulence, and are expected also 
to extend to any desired higher order autocorrelation equations. For example, one may expect 
that the fundamental relations for the third moments in helical turbulence obtained in Chkhetiani 
(1996) will also be expressed simply and concisely in the LANS— a model formulation. This will 
be discussed elsewhere. 

Spectral representation. The spectral representation of these correlation dynamics is also 
worth examining further, especially because the inverse of the Helmholtz operator is algebraic in 
spectral space. As mentioned earlier, this direction was partially developed in Foias, Holm and 
Titi (2001a), where the k'^/^ k'^ roll-off in the LANS —a energy spectrum was discovered in 
the spectral range for which ka < 1 passes to ka > 1. However, much remains to be studied about 
how the alpha- filtering in the LANS— a model affects (or preserves) other fundamental aspects of 
turbulence modeling, while at the same time making the resulting regularized turbulence model 
more computationally accessible. 

Loitsyanski invariant. On the basis of our earlier reasoning that the LANS— a model results 
should recover the NS results as a/r —> 0, we expect that the LANS— a predictions should agree 
with NS theory at large separation, particularly in regard to the behavior of the Loitsyanki 
integral. 

Future directions.^ The KH— a theorem confirms analytically that alpha-filtering leaves the 
velocity autocorrelation statistics invariant at sufficiently large separations a/r <C 1 (and, thus, 
could be expected to preserve intermittency at those separations) for homogeneous isotropic tur- 
bulence. This has three promising implications for future research directions: 

(1) The alpha — > oo limit equation could be used to examine the small-separation effects of alpha 

on the statistics. 

(2) One may study the dependence on the order of the correlation function of the transition from 

the a-dominated to the a-irrelevant regime. This would determine the a-dependence of the 
scaling regions. 

(3) Needing only to compute scales significantly larger than a fixed value of alpha in LANS— a 

should provide more dynamic range for DNS than in the NS case, when alpha vanishes. 
This suggests one might fix alpha in the LANS— a model and use DNS to seek evidence of 
universality at large scales, especially in turbulence decay. 

Acknowledgments. We are grateful to S. Chen, G. Eyink, U. Frisch, R. M. Kerr, J. E. Marsden, 
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